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Abstract 

We study the diffusive scaling limit for a chain of N coupled oscillators. In order to provide 
the system with good ergodic properties, we perturb the Hamiltonian dynamics with random 
flips of velocities, so that the energy is locally conserved. We derive the hydrodynamic equations 
by estimating the relative entropy with respect to the local equilibrium state, modified by a 
correction term. 
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Introduction 

This paper aims at proving the hydrodynamic limit for a Hamiltonian system of N coupled oscilla- 
tors. The ergodic properties of Hamiltonian dynamics are poorly understood, especially when the 
size of the system goes to infinity. That is why we perturb it by an additional conservative mixing 
noise, as it has been proposed for the first time by 011a, Varadhan and Yau [16] in the context of 
gas dynamics, and then in [11] in the context of Hamiltonian lattice dynamics (see e.g. [1], [2], [6], 
[7], [3], [5], [10], [15] for more recent related works). 

We are interested in the macroscopic behavior of this system as N goes to infinity, after rescal- 
ing space and time with the diffusive scaling. The system is considered under periodic boundary 
conditions - more precisely we work on the one-dimensional discrete torus TV := {0, ...,N — 1}. 
The configuration space is denoted by J)jv := (Ix JR) TiV . A typical configuration is given by 
w = {Px, r x ) x< zj N where p x stands for the velocity of the oscillator at site x, and r x represents the 
distance between oscillator x and oscillator x + 1. The deterministic dynamics is described by the 
harmonic Hamiltonian 



N-l r 2 2 
Px I ^jr 

2 2 



(0.1) 

=o 

The stochastic perturbation is added only to the velocities, in such a way that the energy of particles 
is still conserved. Nevertheless, the momentum conservation is no longer valid. The added noise can 
be easily described: each particle independently waits an exponentially distributed time interval and 
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then flips the sign of velocity. The strength of the noise is regulated by the parameter 7 > 0. The 
total deformation ^ r x and the total energy ^2(p x + r x )/2 are the only two conserved quantities. 
Thus, the Gibbs states are parametrized by two potentials, temperature and tension: for /3 > and 
A € R, the equilibrium Gibbs measures /i^ A on the configuration space f2 := (R x ]R) Tjv are given 
by the product measures 

= n - 1 z< d x) drxdpx ' (a2) 

where e x := (p x + r x )/2 is the energy of the particle at site x, and Z(f3,\) is the normalization 
constant. The temperature is equal to (3~ 1 and the tension is given by A//3. 

The goal is to prove that the two empirical profiles associated to the conserved quantities converge 
in the thermodynamic limit ./V — > 00 to the macroscopic profiles r(i, •) and e(t, •) which satisfy an 
autonomous system of coupled parabolic equations. Let ro : T —> R and eo : T — > R be respectively 
the initial macroscopic deformation profile and the initial macroscopic energy profile defined on the 
one-dimensional torus T = [0,1]. We want to show that the functions r(t,q) and e(t, q) defined on 
Mi x T are solutions of 



d t r = -d 2 q v , 

Q te= J_a o 2 f e +- 
2 7 q V 2 



q eT, (0.3) 



with the initial conditions r(0, •) = ro(-) and e(0, •) = eo(-)- 

We approach this problem by using the relative entropy method, introduced for the first time by 
H. T. Yau [19] for a gradient 1 diffusive Ginzburg-Landau dynamics. For non-gradient models, 
Varadhan [18] has proposed an effective approach. Funaki et al. followed his ideas in [12] to extend 
the relative entropy method to some non-gradient processes and introduced the concept of local 
equilibrium state of second order approximation. 

The usual relative entropy method works with two time-dependent probability measures. Let us 
denote by /Iq the Gibbs local equilibrium associated to a deformation profile ro and an energy 
profile eo (see (1.8) for the explicit formula). As we work in the diffusive scaling, we look at the 
state of the process at time tN 2 . We denote it by ^ and we suppose that it starts from . Let 

X e(t,-),r(t,-) 

,N ._ ,,N „ m _ 00 ±N 



^e(t •) r(t ■) ^ e ^ e Gibbs local equilibrium associated to the profiles r(t, •) and e(i, •) which satisfy 
(0.3) 2 . If we denote by and , respectively, the densities 3 of ^ and fi% t n , t ^ with respect 
to a reference equilibrium measure fi^ := Hi , we guess that is a good approximation of the 
unknown density fF. We measure the distance between these two densities by their relative entropy 



H N (t):= I ff(u)log^\d^(u 
Jq n <Pt ( w ) 

Then, the strategy consists in proving that 



lim ^# = , (0.5) 



A conservative system is called gradient if the currents corresponding to the conserved quantities are gradients. 
2 For the sake of readibility, in the following sections we will denote it by ^ t t,) A t ( )i where /?*(■) and At(-) are the 
two potential profiles associated to r(t, ■) and e(t, ■) (see (1.5) and (1.8)). 
3 The existence of these two densities is justified in Section 2.1. 
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and deducing that the hydrodynamic limit holds (for this last step, see [13], [16] or [4]). In the context 
of diffusive systems, the relative entropy method works if the following conditions are satisfied. 



• First, the dynamics has to be ergodic: the only time and space invariant measures for the 
infinite system, with finite local entropy, are given by mixtures of the Gibbs measures in 
infinite volume ug \ (see (1.15)). From [11], we know that the velocity-flip model is ergodic in 
the sense above (see Theorem 1.3 for a precise statement). 

• Next, we need to establish the so-called fluctuation-dissipation equations in the mathematics 
literature. Such equations express the microscopic current of energy (which here is not a 
discrete gradient) as the sum of a discrete gradient and a fluctuating term. More precisely, 
the microscopic current of energy, denoted by j x ,x+i, is defined by the local energy conservation 
law 



where C is the generator of the infinite dynamics. The standard approach consists in proving 
that there exist functions f x and h x such that the following decomposition holds 



Equation (0.7) is called a microscopic fluctuation-dissipation equation. The term Ch x , when in- 
tegrated in time, is a martingale. Roughly speaking, Ch x represents rapid fluctuation, whereas 
V/ x represents dissipation. Gradient models are systems for which h x = with the previous 
notations. In general, these equations are not explicit but we are able to compute them in our 
model (see (A. 16) and (A. 17)). 

• Finally, since we observe the system on a diffusive scale and the system is non- gradient, we need 
second order approximations. If we want to obtain the entropy estimate (0.5) of order o(N), we 
can not work with the measure /i^L x r , t y. we have to correct the Gibbs local equilibrium state 
with a small term. This idea was first introduced in [12] and then used in [17] for interacting 
Ornstein-Uhlenbeck processes, and in [14] for the asymmetric exclusion process in the diffusive 
scaling. However, as far as we know, it is the first time that this is applied for a system with 
several conservation laws. 

Recently, Even et al. [10] used the relative entropy method for a stochastically perturbed Hamilto- 
nian dynamics which is quite close to the dynamics of this paper: the time evolution is governed by 
the same Hamiltonian of anharmonic oscillators but the process is perturbed by a different noise - 
velocities are exchanged and not flipped. Besides, the boundary conditions are mechanical instead 
of periodic. Contrary to this paper, the model is studied in the hyperbolic scale, so that the authors 
do not need to modify the local equilibrium state. 

Up to present, the derivation of hydrodynamic equations for the harmonic oscillators perturbed by 
the velocity-flip noise is not rigorously achieved (see e.g. [7]), mainly because the control of large 
energies has not been considered so far. Indeed, to perform the relative entropy method, we need to 
control the moments 



C'&X — ^ jx— l,x 



(0.6) 



jx,x+i — V/ x + Ch x . 



(0.7) 
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for all k ^ 1, uniformly in time and with respect to N. In fact, the only first several moments are 
necessary to cut-off large energies (as it is explained in Section 2.2) and we need all the others to 
obtain the Taylor expansion which appears in the relative entropy method (see Proposition 2.1 and 
Lemma A. 2). Usually, the entropy inequality (2.24) reduces the control of (0.8) to the estimate of 
the following equilibrium exponential moments 



with 5 > small. Unfortunately, in our model, these integrals are infinite for all k ^ 3 and all 5 > 0. 

To avoid this problem, we could cut-off large velocities by taking a relativistic kinetic energy (as 
done in [16]). Nevertheless, we should change the physics of the problem by modifying the Liouville 
operator, and consequently the fluctuation-dissipation equations would not be available any more. 
Similar difficulties have already appeared in other models: in [8], Bertini et al. do not have these 
precious exponential moments to derive rigorously their results. In an other context, Bonetto et 
al. [9] study the heat conduction in anharmonic crystals with self-consistent reservoirs, and need 
energy bounds to complete their results. Bernardin [2] deals with a harmonic chain perturbed by 
a stochastic noise which is different from ours but has the same motivation: energy is conserved, 
momentum is not. He derives the hydrodynamic limit for a particular value of the intensity of the 
noise. In this case the hydrodynamic equations are simply given by two decoupled heat equations. 
The author highlights that good energy bounds are necessary to extend his work to other values of 
the noise intensity. In fact, only the following weak form is proved in his paper: 



In this work, we get uniform control of (0.8) (Theorem 1.2), thanks to a remarkable property of our 
model: the set of convex combinations of Gaussian measures is preserved by the dynamics. This is 
one of the main technical novelties in our work. 

The next section contains a more precise description of the results outlined here, along with the plan 
of the paper. 

1 The Model and the Main Results 
1.1 Velocity-flip Model 

We consider the unpinned harmonic chain perturbed by the momentum-flip noise. Each particle has 
the same mass that we set equal to 1. The configuration space is denoted by Q N := (R x M) Tjv . 

A typical configuration is u = (r,p) € £l N , where r = {r x ) X £j N and p = (p x )xeJ N - 

The generator of the dynamics is given by £n := An + ^Sn, where, for any continuously differen- 
tiable function / : Q N — > M, 





(0.9) 



Ar(/) := [(Pz+l - Px) d r J + (r x 



r x -i) dp J] 



(1.1) 
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and 



S N (f)(r,p) := \ £[/(r,p*)-/(r,p)] 



(1.2) 



Here p x is the configuration obtained from p by the flip of p x into — p x . The parameter 7 > 
regulates the strength of the random flip of momenta. 

The operator An is the Liouville operator of a chain of interacting harmonic oscillators, and Sn is 
the generator of the stochastic part of the dynamics that flips at random time the velocity of one 
particle. The dynamics conserves two quantities: the total deformation of the lattice TZ = YlxeJ N r x 
and the total energy £ = YlxeT N e x , where e x = (p x + r 2 )/2. Observe that the total momentum is 
no longer conserved. 

The deformation and the energy define a family of invariant measures depending on two parameters. 
For (3 > and A € R, we denote by fJ^ x the Gaussian product measure on £l N given by 



,—Pe x —Xr x 



-dr x dp x . 



An easy computation gives that the partition function satisfies 

2n 



2p) 



(1.3) 



(1.4) 



In the following, we shall denote by //[•] the expectation with respect to the measure fj,. We introduce 
L, 2 (pp A ), the space of functions / defined on f2 such that A^aI/ 2 ] < +00. This is a Hilbert space, 
on which An is antisymmetric and Sn is symmetric. 

The thermodynamic relations between the averages of the conserved quantities f E R and e € 
(0, +00), and the potentials (3 E (0, +00) and A £ R are given by 

A 2 



e(/3,A) :=$ tX [ 



f(/3,A) :=^ x [r x 



P + 2/3 2 
_A 

~1 ' 



Let us notice that 



V (3 e (0,+oo),V A GR, e(/3,A) > 



r 2 (/3,A) 



(1.5) 



(1.6) 



We assume that the system is initially close to a local equilibrium (defined as the following). 

Definition 1. A sequence (^i n )n of probability measures on Q N is a local equilibrium associated 
to a deformation profile tq : T — > R and an energy profile eo : T — > (0, +00) if for every continuous 



function G : T 



and for every S > 0, we have 



lim u 



N 



lim fi 

iV->-oo 



N 



xGTm 



r x 



e x 



G(q)r (q)dq 
G(q)e (q)dq 



> 5 



> 5 



, 



. 



(1.7) 
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Example 1.1. For any integer N we define the probability measures 

u N (dr dv) - TT e M-Po(x/N)e x - X (x/N)r x ) 

where (3q(-) an d Ao(-) are related to eo(-) and ro(-) by (1.5) 

= §(#,(■), A (-)) , 




= r(A)(-),A (-)) • 

Then, the sequence (/^q a ( )) * s a ^ occ ^ equilibrium, and it is called the Gibbs local equilibrium 
state associated to the macroscopic profiles eo and tq. Both profiles are assumed to be continuous. 

To establish the hydrodynamic limit corresponding to the two conservation laws, we look at the 
process with generator N 2 £jy, namely in the diffusive scale. The configuration at time tN 2 is 
denoted by uj^ , and the law of the process (uj^)t^o is denoted by nf . 



1.2 The Thermodynamic Entropy 

The function 

5(e,r)= inf {Ar + /3e + log Z(f3, A)} (1.9) 

/3>0,AGM 

is called the thermodynamic entropy. An easy computation, coming from the explicit expression of 
the partition function, gives 

/ r 2 \ r 2 
5(e, r) = 1 + log(27r) + log I e - — J , when e - — > . (1.10) 

The relations (1.5) can be inverted according to 

dS(e,r) dS(e,r) 
(e ' r) = dr ' ^ (e,r)= de " ( ^ 

Remark. These two equalities, together with (1.5), show that there exists a bijection between 
the two sets {(/3,A) G M 2 ; /3 > 0} and {(e, r) £ I 2 ; e > r 2 /2}. From the equations above, the 
inverted relations can be written as 

A(e,r) = -^-^, PM = -±-. (1.12) 

We denote by \E' the function 

*: {(e,r) G M 2 ; e > r 2 /2} -> {(/3, A) G M 2 ; /3 > 0} 
(e, r) ^ 



ie -r 2 /2' e-r 2 /2j ' 
If ?7 = (e, r) and x = (A A) satisfy the relations (1.5), then n and x are sa id duality and we have 

-5(e,r) + logZ(^A) = -7 ? - X - (1.13) 
Here, the notation a ■ b stands for the usual scalar product between a and b. 
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1.3 Hydrodynamic Equations 



Let [i and v be two probability measures on the same measurable space (X, J 7 ). We define the 
relative entropy H([x\v) of the probability measure fi with respect to the probability measure v by 



H{n\v) = sup / dji - log e f di^J j 



;i.!4) 



where the supremum is carried over all bounded measurable functions / on X. 

The Gibbs states in infinite volume are the probability measures fis,X on f2 = (M. x M) z given by 

-j3e x -\r x 



(dr,dp) = ] 



xG2 



Z(/3,A) 



-dr x dp x 



;i.i5) 



We denote by r x c/? the shift of (t x (/?)(u;) = (/9(r x u;) = + • )). In this article the following 

theorem is proved. 

Theorem 1.1. Let (^q)n be a sequence of probability measures on Q N which is a local equilibrium 
associated to a deformation profile tq and an energy profile eo such that eo > Iq/2 (see (1.8)). We 
denote by /3q and Aq the potential profiles associated to tq and eg.' 



(A),A ) := *(e ,r ) 



We assume that 



1.16) 



and that the initial profiles are continuous. 

We also assume that the energy moments are bounded: let us suppose that there exists a positive 
constant C which does not depend on N and t, such that 



v k ^ 1, 



N 



E4 



[Ckf N . 



;i.i7) 



Let G be a continuous function on the torus T and if be a local function which satisfies the following 
property: there exists a finite subset A C Z and a constant C > such that, for all oj G Vt N , 
<p(u) < C (1 + Ei G A e i( w ))- Then > 



N 



Af-s-oo 



-> 



(1.18) 



where (p is the grand- canonical expectation of (p: in other words, for any (e, r) € M 2 suc/i i/ioi 
e > r 2 /2, Zei (/3, A) = #(e,r) i/ien 



J(RxR) z 



;i.i9) 
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Besides, e and r are defined on M + x T and are solutions of 

1 

"f 



d t r = -dgT , 



dte = ±% | 



g t, t g 



;i.20) 



o(iV). 



with the initial conditions r(-,0) = ro(-) and e(-,0) = eo(-)- 
Remarks. 

1. In order to prove the theorem, we shall show afterwards that H (t^t^l^^.) 

Here ^w.) is a probability measure which is close to the Gibbs local equilibrium fi^ t ^ w t n 
(1.8). The functions (/3(t, •), A(t, •)) are still related to e(t, •) and r(t, •) by (1.5). 

This fact allows to establish the hydrodynamic limit in the sense given in the theorem. For a 
proof, we refer the reader to [16], Corollary 2.2, [13] or [4]. 

2. Let us notice that the functions e, r, (3 and A are smooth when t > 0, since the system of 
partial differential equations is parabolic. Moreover, the function = e — r 2 /2 satisfies 



1 



0, 



1.21 



The supersolutions of the heat equation follow the minimum principle. Consequently, since 
there exists c > such that the initial profile /?o has the following property 

V q G TV, (3 (q) ^ c> , 

then we know that the function j3 satisfies: 

VqeT N , Vte[0,T], /5t(g)>c>0. (1.22) 

3. After some integrations by parts, a simple computation shows that 

+ 2p(t,q) [d q r(t,q)] 2 dq^O (1.23) 



} 27/ 









when r and e are the solutions of the hydrodynamic equations (1.20). This fact is in agreement 
with the second thermodynamic principle. 

In Section 3, we will show that the hypothesis on moments bounds (1.17) holds for a wide class of 
initial local equilibrium states. Before stating the theorem, let us give some definitions. 

We denote by ©at(R) the set of real symmetric matrices of size N. The correlation matrix C G 
S2At(K) of a probability measure v on Q N is the symmetric matrix C = {Ci j)l<gi,j<2iV defined by 

v[nrj] i,j G {1, ...,N} , 

u[r iPj ] i G {1, N}, j G {N + 1, 2iV} , 

u\p irj ] i G {N + 1, 2N}, j G {1, N} , 

vfriPj] i,j e{N + l,...,2N] . 



a 



1,] •" 



(1.24) 
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Let us denote by E^r the subset of M? N x (5 



2N\ 



defined by the following condition: 



(m, C) G Sat <^> < 



m k 
Cu > 





= 



c. 



7717 



i+N,i+N 



for all it = iV + 1 . . . 2N 
for all i j , 
for all i = 1 . . . 2N , 
for all i = 1 . . . N . 



(1.25) 



Precisely, it means that m is of the form m = (mi, . . . , mjv, 0, . . . , 0), and C is a diagonal matrix 
whose diagonal components can be written as (mf + a\, . . . , + oat, ai, . . . , ayv), where cti > 
for all i = 1 . . . N. 

For (m, C) € S^r, we denote by G m ,c(') the Gaussian measure with mean m and correlations given 
by the matrix C. The covariance matrix of G m fi{') is thus C — mtm. 

Lemma 1.1. Let A(-) anc? /3(-) 6e too functions of class C 1 defined on T, and be the Gibbs 

local equilibrium defined by (1.8). If we denote by mpr.yu.) an d C/3(-),A(-) respectively the mean vector 
and the correlation matrix of A(-)' then we have 

(m ( g ( . )iA( . ) ,C / g ( . )>A{ . ) ) G Sat and M/»(.),A(-) = G ^ ( .),a ( .),C / , ( . )iA(0 - 

Proof. This result comes from the explicit formula of ^go A o given in (1.8). First, notice that each 
momentum p x is centered under _w \ and 



/^(.),A(.)M " (]v) 
Second, we easily obtain the following expressions: 



m /3('),A(-) 



A f0_ 

1 \N 




(1.26) 



[1.27) 



O^v D' 



where < 



D = diag 
D' = diag 



1 ^ A 2 (x/iV) 



' P(x/N) p 2 {x/NY 
1 



/3(x/iV) 



(1.28) 



□ 



Now we state our second main theorem, which will be proved in Section 3. 

Theorem 1.2. We assume that the initial probability measure Hq is a convex combination of Gibbs 
local equilibrium states. More precisely, let a be a probability measure whose support is included in 
Sat. We assume that a satisfies: 



for all k 



^ 1, J [K(m,C)] h da(m,C) < oo , 



(1.29) 
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where K(m,C) := sup {Co}. We define the initial probability measure fj,^ by 
i=l...N 

= J Gm,c{-) da(m,C) . (1.30) 
Then, (1.17) holds, and the conclusions of Theorem 1.1 are valid. 

Remark. As in [7], we could consider a more general model, with a pinning potential. Instead of 
the deformation r x , we now introduce the position q x of the particle x. The new pinning Hamiltonian 
is given by 

«5r=Ef+" 2 £f+ E ^T^- (1-31) 

x-eTjv rreTjv \x—y\=l, 

x,y€t N 

The strength of the pinning potential is regulated by the parameter v > 0. The energy of site x is 
now given by 



l+^f+i E (1-32) 



2 ' 4 

y;\x-y\=l 



The stochastic operator remains equal to Sjv, and the Liouville operator J? N can be written as 
follows: 

A% = E ^ - ^ - > (1-33) 

where A is the discrete Laplacian: (Au) x = u x +x + u x ._i — 2u x . 

Because of the presence of the pinning, the bulk dynamics conserves only one quantity: the total 
energy e x . It follows that the Gibbs equilibrium measures pJS are fully characterized by the 
temperature Under fi^, the variables p x are independent of the q x and are independent 

identically Gaussian variables of variance /3" 1 . The q x are distributed according to a centered 
Gaussian process with covariances given by 

$ (q x q y ) = T{x - y), such that [(u 2 - A)] T(z) = -l z=0 • (1-34) 



Observe that there exists C := C(v) independent of N such that /J,g (q x q y ] 
N>1. 



^ C- l e~ c \ x - y \ for 



any 



These correlations make computations more technical, but the hydrodynamic limit can be established 
by following the proof here (in [7], Section 3.2, a heuristic argument is given). Assume that the 
system is initially distributed according to a Gibbs local equilibrium associated to the energy profile 
eo(-)j and define e(t, •) as the evolved profile in the diffusive scale. Then, if the energy moments are 
bounded like (1.17), e is the solution of the following heat equation 

{dte = d„(D(e)d a e) , 
q{ ; q ' qeT, tett, (1-35) 
e(0,-) = e (-) , 

where -D(e) is the diffusivity given by 

D := Die) = ^ 7 . (1.36) 
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In our model, where the state space is not compact, what matters is the existence of moments 
bounds. We will see in Section 3 that this existence can be easily justified by following the same 
ideas which work for the unpinned model. 

For the sake of simplicity, we will denote by et(-), rt(-), \t(-) and Pt(-), respectively, the functions 
q — > e(t, q), q — > r(i, q), q — > X(t, q), and q — > f3(t, q) defined on T. 

1.4 Ergodicity of the Infinite Velocity-flip Model 

We conclude this part by giving the theorem of ergodicity, which is proved in [4], Sections 2.2 
and 2.4.2, by following the ideas of [11]. Let us define, for all finite subsets A C Z, and for two 
probability measures v and /i on f! := (R x the restricted relative entropy 

H A {v\n) := H{u A M (1.37) 

where v\ and //a are the marginal distributions of v and /ion ft. 

The Gibbs states in infinite volume are the probability measures /x^ \ on ft given by 

e ~Pe x -\r x 

d w,>> ; = n z(b \) drxdpx • ^- 38 ^ 

The formal generator of the infinite dynamics is denoted by C 

Theorem 1.3. Let u be a probability measure on the configuration space ft such that 

1. v has finite density entropy: there exists C > 0, such that for all finite subsets AcZ, 

H A (u\[j,*) ^ C\A\ , (1.39) 
with /x* := /ii 5 o a reference Gibbs measure on (M x M) z , 

2. v is translation invariant, 

3. v is stationary, i.e. for any compactly supported and differentiable function F(r,p), 

J A(F) dv = . (1.40) 

4- the conditional probability distribution of p given the probability distribution ofr, denoted by 
u(p\r), is invariant by any flip p — > p x , with x € Z. 

Then, v is a mixture of infinite Gibbs states. 

Corollary 1.1. If v is a probability measure on ft satisfying 1, 2 and if u is stationary in the sense 
that: for any compactly supported and differentiable function F(r, p), 

J C{F) dv = , (1.41) 

then v is a mixture of infinite Gibbs states. 
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The outline of the rest of the paper is as follows. In the next section we expose the strategy of the 
proof. We introduce the relative entropy -ffjv(i) of /if with respect to a corrected local equilibrium, 
and we prove a Gronwall estimate of the entropy production of the form 

d t H N (t) ^ C H N (t) + o(N) , (1.42) 

where C > does not depend on N. In Section 3 we prove Theorem 1.2. 

We suppose that t belongs to a compact set [0, T], T fixed. All estimates are uniform in t G [0, T]. 



2 Entropy Production 

2.1 Introduction to the Method 

For the sake of simplificity, we denote all couples of the form (/?(•), A(-)) by \{' 
The corrected Gibbs local equilibrium state v^ t /.\ is defined by 



^« := zmy) n ex p GO e * - a * (j?) r * + ^ f (*' |) • r - h ^ p)) dr * d p* 



where Z(xt(')) is the partition function and F, h are functions which will be precised later on. The 
notation a ■ b still stands for the usual scalar product between a and b. An estimate of the partition 
function Z(xt(')) is performed in Appendix A. 

We are going to use the relative entropy method, with the corrected local Gibbs state i.\ instead 



of the usual one /j, n , v We define 



xt(-) 

*(•)■ 



H N (t) := H („f |^ ( .)) = j nN ft N («)log^g^ ( W ) , (2.2) 

where is the density of /if with respect to the reference measure /if - This is a solution, in the 
sense of the distributions, of the Fokker-Planck equation 

d t f t = N 2 C* N f t (2.3) 



where C* N = —An + 7SV is the adjoint of Cn i n L 2 (/if ). In the same way, is the density of 

xt(0 w ^ n res P ec t t° (which here is easily computable). 
Thus, our purpose is now to prove (1.42). We begin with the following lemma. 



Lemma 2.1. 

r 1 ... .. r 1 „ ... 1 

(2.4) 



d t H N (t) ^ J -L {N 2 C* N ^ - dttf) f t N d^ lfi = /if 



Proof. The case where is smooth is proved in [13], Chapter 6, Lemma 1.4. Here, we do not 
know that is smooth, so that we refer the reader to the proof in [5], Section 3.2, which can be 
easily followed. □ 
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Now, we choose the correction term. We consider 

x 

< 



T x h(r,p) :-- 



N 



27 



Px+1 + Px + 7r r a 



Px+1 



(2.5) 



Thus, 



f (r,p) 



(z(i,o)) n 
z(xt(-)) 



n 



AT 



x 



exp -A — + 1 - A t — )r x + —F [t, — ■ r^r, p 



A r 



1 



N 



N 



(2.6) 

We define £ x := (e x ,r x ) and r](t,q) := (e(t, q),r(t, q)). If / is a vectorial function, we denote its 
differential by Df. 



In Appendix A, the following technical result is proved. 

if)- 1 (AT^f 

a microscopic expansion up to the first order appears. In other words, 



Proposition 2.1. The term (4>t) 1 (yN 2 C* N (f)f — dt<ft^) is given by the sum of five terms in which 



1 



(AT^f - d t J>?) 



fc=l xSTjv 



4~H k V t 



N 



(DH k )( V (t, 



N 



N 



+ o(N) (2.7) 



where 



k 


J x 


H k (e,r) 


Vk(t,q) 


1 


pi + + 2^p x r x _i 


r 2 

e H 

2 


^d 2 q P(t,q) 

--d 2 X(t,q) 

7 


2 


r x + jp x 


r 


3 


pl (r x + r x -i) 2 


(2e-r 2 ) (e + V) 


^[W,g)] 2 


4 
5 


P 2 X (r x + r^-i) 

pI 


r (2e - r 2 ) 
r 2 

e 

2 


-d q /3(t,q) d q X(t,q) 
7 

^[^A(i,g)] 2 



A priori the first term on the right-hand side of (2.7) is of order N, but we want to take advantage 
of these microscopic Taylor expansions. First, we need to cut-off large energies in order to work 
with bounded variables only. Second, the strategy consists in performing a one-block estimate: we 
replace the empirical truncated current, which is averaged over a microscopic box centered at x, by 
its mean with respect to a Gibbs measure with the parameters corresponding to the microscopic 
averaged profiles. 

A one-block estimate will be performed for each term of the form 



X 

N 



x 
N 



(2.8) 
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In the following the index k is omitted, whenever this does not cause confusion. We follow the lines 
of the proof given in [4], Section 3.3 and inspired from [16]. A sketch of the proof for the one-block 
estimate is given in Appendix B. 



2.2 Cut-off of Large Energies 

For x G T N , we define A X)M ■= {e x + e x -\ ^ M}, J XjM ■= Jx ^A XiM , and £ X)M := £ x t e ^M- 

Then, J X ,M an d £ X ,M are bounded by C{M) > 0. 

We use twice the Cauchy-Schwartz inequality to write 

\ 1/2 / \ 1/2' 

1/2 



JV 



E» 2 (^ 



/'/ 



N 



y i A ° 



1/2 



(2.9) 



First, v 2 (t,x/N) is bounded by a constant which does not depend on A. Second, the term J 2 can 
be bounded above by the squared energy e 2 . The hypothesis (1.17) shows that there exists Co which 
does not depend on A such that 



£" 2 (^ 



1/2 



C N l > 2 . 



(2.10) 



Moreover, Markov inequality proves that 



Finally, we obtain a constant C independent of A such that 



N 



< CA e(M) . 



(2.12) 



Observe that this estimate is in agreement with the Gronwall inequality we want to prove, since we 
are going to divide by A. Thus, the error term is of order 1/M that goes to as M — > oo. 

Consequently, J x can be replaced by J X ,M i n (2-8), and similarly, £ x can be replaced by £ X ,M- 
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2.3 One-block Estimate 



Now we prove that 



N 



£ v (t, I) [J.^ - H („ (t, I)) - (DH) (p (t, £)) • fa - „ (t, £))] 

H N (t) 



<C -j^ + £ (N,M) (2.13) 



with e(N, M) 

We dem 
profiles 



-> 0. 



A/— >oo,A r — >oo 

We denote by A^(y) the box of length £ centered around y. We introduce the microscopic average 



3'GA £ (y) 



(2.14) 



We split Ttv into p = iV/-£ boxes A^(xj) centered at Xj. Here £ is assumed to divide N for simplicity. 
We will first let A" — > oo, then £ — >• oo and then M — >• oo. 



First of all, we want to replace 



by 



3=1 



7 £ J ^ 

iGA^(xj) 



(2.15) 



(2.16) 



The error term produced during this step can be written as 



E E 

3=1 ieA^ay) 



J, 



i.M 



N J V N 

The last inequality comes from the smoothness of v, more precisely 



< Cx(M)- . 



(2.17) 



(2.18) 



Similarly, we perform the same estimates for the other terms and it remains to prove that 



P 



3=1 



at; £ 



A 



N 



Ve,M( x j) -V (t, 



(2.19) 



vanishes as M, A", £ — >■ oo, the limit in A" taken first, then the limit in £ and finally the limit in M. 
The additive term which appears after performing this replacement can be bounded above by a term 
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e jv m,i which depends on N, M and £, but which is independent of the particular splitting of HV 
into p boxes. This term is of order o(N) in the Gronwall inequality we want to prove, in the sense 
that 

. (2.20) 



lim lim lim N 1 fif [e n m i] 

M— >oo I— yoo N— >oo 



Now we want to perform a one-block estimate. The main idea consists in replacing £ 1 Yli<=A e (x ) ^ 



•i.M 



by H(rj£^M(xj)). This is achieved thanks to the ergodicity of the dynamics (see Theorem 1.3). In 
order to use this ergodicity property, we have to work with a space translation invariant measure. 
To obtain such a probability measure, we introduce a second average over the Xj, 1 ^ J ^ p. For 
each k € {0, £ — 1}, we can split Tjv into p disjoint boxes of length £ by writing 



V k = 0, ...,£- 1, T N = \J At(xj + k) 



(2.21) 



Then, we average the different splittings mentioned above. More precisely, in Appendix B we recall 
how to prove 



1 l ~ X 

lim sup lim sup lim sup - ji^ 



M— »oo I— >oo JV— »oo 



fc=0 



1 P 



v t 



[xj + k] 
N 



' ir- A .f™ . I U\ 



ieA e (xj+k) 



= 0. 
(2.22) 



2.4 Large Deviations 



The previous estimates are valid for any splitting of Tjy into p boxes of length £. Thus, it would be 
sufficient to prove (2.19) with every Xi replaced by Xi + k for arbitrary k £ {1, ...,£ — 1}. Consequently, 
it is sufficient to prove (2.19) in an averaged form. Then, from the one-block estimate, we have to 
deal with 



1 



k=0 



N 



1 



[ x j + k \ 
N 



[ x j + k ] 
N 



(2.23) 



where 0(w, u) := H(w) - H(u) - DH(u) • (w - u). 

By definition of the entropy, for any a > and any positive measurable function / we have 



J fdfi^^ {log (J e a f du^j + H(n\v)\ 



(2.24) 



This inequality, known as the entropy inequality, allows to show that: for any a > 0, (2.23) is less 
than or equal to 



HN(t) 1 1 , N 



a N 



k=0 



(2.25) 



Notice that the last integral converges because all quantities are bounded. 
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The first term is in agreement with the Gronwall inequality we want to obtain. We look at the second 
term. Since we have arranged the sum over p disjoint blocks which are independently distributed 
by v ^ t (.y the second term is equal to 



1 P 

k=0 3=1 



<(■) 



nt,M{xj+k),r)[ t 



[xj+k] 



(2.26) 



We are going to show that this expression vanishes as M, N, £ — > oo by using the large deviation 

,,N 



properties of the measure u N ,s, that locally is almost homogeneous. In fact, by using the smoothness 



for the various involved functions, we can substitute the inhomogeneous product measure u 



»(•) 



restricted to Ag(xj + k) with the homogeneous product measure M^ t ([x ■+&]/#)> in eacn expectation 
of the previous expression. More precisely, we have the following lemma. 

Lemma 2.2. 



M x {N,l,k,M) : =-L^log 
can be replaced by 



af 



t> 1 1, ^ '" J I n I m ,m (xj+k) ,r) 1 1, [ ^ fcl J 



(2.27) 



1 p 

M 2 (N,£,k,M) : =_]Tlog^ 



*»(fo+k]/tf) 



j'=i 



of 



;,■+*] 



)) 



(2.28) 



TTje difference between these two terms is less than or equal to a small term which depends on I (but 
not on k) and vanishes in the N limit: there exists a constant C(£,M,N) which does not depend on 
k such that 

MAN, £, k, M) - M 2 (N, £, k, M) ^ C(£, M, N) and C(£, M, N) ► . (2.29) 

N— >oo 



Remark. In the following, we will prove that 



lim sup lim sup lim sup M 2 (N, £, k, M) = . 

M— >oo I— >oo N— >oo 



In addition to this lemma, this implies that 



lim sup lim sup lim sup Mi (N, £, k, M) = , 

M— >oo ^— >oo N-^oo 



(2.30) 



(2.31) 



since M±(N,£,k,M) is always nonnegative, and we know that, for all sequences a n and b n , 

lim sup a n ;C limsup(a n — b n ) + lim sup b n . (2.32) 
Proof. For each j € {1, ...,p}, the function 



Fj := exp < a£ 



vlt, 



[Xj + k] 
N 



^ Vi,m( x j + k),rj(t, 



[xj + fc] 
N 



(2.33) 
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is bounded above by e , C > (since r]£M i s bounded and t belongs to a compact set), and 
depends on the configuration only through the coordinates in Ap(xj + k). Thus, each expectation 
appearing in the sum can be taken w.r.t the restriction to Ag(xj + k) of u ^ t r.y These restrictions 
are inhomogeneous product measures but with slowly varying parameters and hence, each term 



log v^ t ^[Fj] can be replaced by log// 



N 



[FA with a small error. 



Indeed, the difference between these two terms is equal to 



i N 

lo ZVx t (\x 



Xt([xj+k]/N) 



1 + 



Fj (hj - 1) 

V Xt ([ X j+k]/N)i F A 



^Xt([x 3 +k]/N)i F j (hj 1)] 



log 1 + 



N 



Xt(lxj+k]/N)l F l 



(2.34) 



with 



hj := exp I ^2 €i,M ■ 

v i£A e (xj+k) 



Xt 



N 



Xt 



X j ~\~ k 

N 



N 



Fit, 



N 



nh + 



log Z xt 



+ 
i 

N 



log Z xt 



N 



(2.35) 



The inequality log(l + x) ^ \x\ (true for any real x) and the fact that M^([ x +k]/N) l x 3 



us to estimate 



Vxt{[x 3 +k\/N)[\ F i( h i ~ !)|] 



[Fj] ^ 1 reduces 
(2.36) 



By using the smoothness of Xt an d the inequality | e 33 — 1 1 ^ |x|e' x ', one easily shows that there exist 
positive constants Co, C(£), and /3 which do not depend on j such that 



1^-1)1 < 



C(£)£ 
N 



^ [e-i,M + e i+lsM + 1] J exp 



Co£ 
N 



^ [ei,M + ej+i,M + 1] f , 

i€A e (xj+k) 



and 



d ^Xt{[xj+k}/N) 



< C(i) . 



(2.37) 



Ai(xj+k) 



Hence, the total error performing by these replacements is bounded above: 



exp < 



ieAf(o) 



for some positive constant Ci(£,M). 

It trivially goes to as N goes to infinity for each given fixed I. 



□ 



Lastly, we have to show that the limit 

i-l -, p 



11 

lim sup lim sup lim sup - ]T — £ log H$ t{[xj+k]/N) 

fe=0 j=l 

vanishes. Here, rji iM ■= V£,m(°) = ^ E; 6 A £ (o) &,Af ■ 



[a;,-+fc] 



,-+fc] 



(2.38) 
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The limit in p results in an integral over T because we have a Riemann sum. Moreover, the integral 
does not depend on k so that the averaging over k disappears in the p limit. Hence, the point is to 
estimate 

limsuplimsup^ / log M \ e a£v ^ n( m , M Mt,<i))] dq . ( 2 . 39 ) 
A/^oo i^oo al J T *«w L J 

According to Laplace- Varadhan theorem applied to these product measures an d according 

to the dominated convergence theorem, the previous limit is equal to 

limsup- / sup {av(t,q) n(z,7](t, q)) - I M (z,ri(t,q))}dq , (2.40) 

A/->oo a JT zGK 2 

where I][i(z,rj(t,q)) is the rate function of the sequence |/c _1 X^i=i £«,A/| as ( r x,Px)xeT N are dis- 
tributed according to the homogeneous product measure • 

The function Im is the Legendre transform of the cumulant-generating function of £o,Af : 

I M (z, V (t,q)) = sup {y • z - log ^ t{q) [e y<0 ' M ]} ■ (2.41) 

Hence 

liminfI M (z, ??(*,<?)) ^ sup {yz-log/A Je^ }} = I(z lV (t,q)) , (2.42) 

where I(z,rj(t,q)) is the rate function of |^ _1 X^i=i£«} as { r yiPy)y are distributed according to 
the homogeneous product measure • 

It follows, by Fatou's lemma, that (2.40) is smaller than or equal to 

- / sup {av(t, q) n(z, rj(t, q)) - I(z, rj(t, q))} dq . (2.43) 
a J j z 

From now on we omit the dependance in (t, q) of the involved functions v and t]. Recall that \ an d 
r\ are in duality (see (1.13)). An easy computation gives that 

I{z, rj) = sup jy • z - log (J e y< e x< - l ° sZ(x) drdp^j X 

= sup{y • z - log Z(x + y) + log Z(x)} 
y 

= bgZ( x ) + z- X -S(z) , (2.44) 

where the last equality follows from the equality between the Fenchel- Legendre transform of log Z 
and the function —5. We observe that 1(77,77) = and D z I(z,rj) = 0. Furthermore, / is strictly 
convex in z: 

(Dll)(z, V ) = (Dl{-S})(z)>0. (2.45) 
Since ^(77,7?) = and (D z O)(z,r?) = {DH){z) - DH{rf), we also get: (D z n)(rj, rj) = . 
Lemma 2.3. For a > sufficiently small, 

Vz£K 2 , V q € T, av(t, q)fl(z, i](t, q)) < J(z, rj(t, q)) . (2.46) 
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Proof. An easy computation provides an explicit expression for the rate function: if z = (zi, Z2) and 
77 = (e, r) with e — r 2 /2 > then 

From the inequality — logs ^ — x + 1 (satisfied for any x > 0). we get 

> (r - 22)2 ■ (2 - 48) 

Thus, for a given 77, the rate function z — > J(z, 77) is such that /(z, 77) ^ c^|z — r/| 2 , where is a 
positive constant. Moreover, according to (1.22), 

VtG[0,3l ) V?€T ) c, M >c>0 (2.49) 

Let us fix z € R 2 . From the Taylor-Lagrange theorem, there exists a positive constant C such that 

fi(z, »/(<,<?)) ^C|z-r/(^)| 2 <I(z, »/(*,</)) . (2.50) 

More precisely, C is equal to 

sup \\D 2 H( V (t,q))f . (2.51) 
(t, 9 )e[0,T]xT 

Since v is uniformly bounded, the result is proved. □ 

Consequently, for a small enough, 

sup {av(t, q) n(z, g)) - /(z, r?(t, g))} = , (2.52) 

z 

and we have finally proved that 

d t H N (t) < C ffiv(t) + R N ,l,M{t) (2.53) 

with 

1 /"* 

lim lim lim — / R N 1 M (s) ds = . (2.54) 
By Gronwall's inequality we obtain: H^{t)/N > and Theorem 1.1 is proved. 

N-*oc 

3 Proof of Theorem 1.2: Moments Bounds 

In the following, we prove the two conditions on the moments bounds for a class of local equilibrium 
states. First, we assume that the initial law /j,q is exactly the Gibbs local equilibrium measure 
^A>(-) A (-)' Second, we extend the proof to the case where /Jq is a convex combination of Gibbs local 
equilibrium measures. 

We need to control the moments [^2 X e^] for all k ^ 1. The first two bounds (A; = 1, 2) would be 
sufficient to justify the cut-off of the currents, but here we need more because of Lemma A. 2 (which 
is necessary to prove Proposition 2.1). Since the chain is harmonic, Gibbs states are Gaussian. We 
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recall that all Gaussian moments can be expressed in terms of variances and covariances. In the 
following, we first give an other representation of the dynamics of the process, and then we prove 
the bounds and precise their dependence on k. 

Let us highlight that, from now on, we consider the process with generator Cn' it is not accelerated 
any more. The law of this new process (Cot)t^o is denoted by pf . At the end of this part, Theorem 
1.2 will be easily deduced since all estimates will not depend on t, and the following equality still 
holds: 

^ = ■ (3.1) 

Remarks. 



1. In the following, we always respect the decomposition of the space f2 



N 



R N . Let us 



recall that the first N components stand for r and the last N components stand for p. All 
vectors and matrices are written according to this decomposition. Let v be a measure on Q N . 
We denote by m G R 2N its mean vector and by C G 5D?27v(^) ^ s correlation matrix (see 
(1.24)). We can write m and C as 



m = (p, vr) G R 2N and C = (Z y) 



where p := u[r] G U. N , ir := u[p] G R N and U, V, Z G Wl N (R) . 
2. Thanks to the convexity inequality (a + b) k ^ 2 k ~ l (a k + b k ), for a, b > 0, we can write 



(3.2) 



Thus, instead of proving (1.17) we will show 



ft 



N 



P. 



2 k 



< (Ckf N and p,\ 



N 



E 



,2k 



s: {Ckf N . 



(3.3) 



(3.4) 



3.1 Poisson Process and Gaussian Measures 



We are going to use a graphical representation of the process (uJt)t^o- 
Let us define 



A 



/0 







V(°) 







1 

-1 ■■. 



-1 1 



(0)\ 







: '•• 1 

10 0-1 

-1 

: : 



-1 1 



G Wl 2N {R) 



(3.5) 
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We now consider (mj, Ct)t^o, a Markov process on M? N x ©2jv(K) whose generator is denoted by Q 
and defined as follows. 

(U Z*\ 
2 y j G ©2Ar(lR), where p, it are two vectors in M. N , U, V are 

two symmetric matrices in ©at(1R) and Z is a matrix in S0tjv(]R). Hereafter, we denote by Z* the 
transpose of the matrix Z. 

The generator Qn is given by 



(g N v)(m, C) := (IC N v)(m, C) + 7 (H N v)(m, C) , (3.6) 

k n ■■= (~ AC + CA ^ dc ^ + 5^ - ^) d Pi + - ft-i) 5 ^} > ( 3 - 7 ) 

(H N v)(m, C):=^J2 M m *' C ^ ~ v( - m > C ^ ' ^ 

k£t N 

m k = (p, vr fc ) and C k = • C • £ fc = . (3.9) 

In these last two formulas, ir k is the vector obtained from 7r by the flip of 7Tfc into — tt^, and is 
defined as 

In n 



where 



and 



Here. 



_0 n I n — 2Ek,k J 
More precisely, 

Z k j = {-if^Zij and V k j = (-l)^.i+ s ^)v i j . (3.11) 

We denote by P m o,Co the law of the process (m t , Ct)t^o starting from (mo, Co), and by E m0i c [-] the 
expectation with respect to P m0i c . 

For t ^ fixed, let 6^ q (•, •) be the law of the random variable (mj, Ct) G M. 2N x ©2At(M), knowing 
that the process starts from (mo, Co). 

Recall that we denote by G m c(-) the Gaussian measure on Q N with mean m G M. 2N and correlation 
matrix C G 6 2 at(P)- 

Lemma 3.1. Let Pq := pp Q ^ x (-) ^ e ^ e Gibbs equilibrium state defined by (1.8), where Xq(-) and 
Po(-) are the two macroscopic potential profiles. 

Then, 

ftt= J G m ,c(') <,c„KC) (3.12) 
- ! -|'-^(^),".-|(^).0-.q| (-13) 
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and 



Proof. We begin with the graphical representation of the process (u>t)t^o, which is based on the Harris 
description. Let (Ni)i^j N be a sequence of independent standard Poisson processes of intensity 7. 
In other words, we put on each site i G T/v an exponential clock of mean I/7. At time the process 
has an initial state uq. Let T\ = inf^o {3 i G HV, 2Vj(i) = 1} and i\ the site where the infimum is 
achieved. 

During the interval [0, Ti), the process follows the deterministic evolution given by the generator 
An- More precisely, let F : (r, p) G — > A - (r, p) G where A is given by (3.5). Then, for any 
continuously differentiable function / : Q N — > R, 

A N f{oj) = A ■ Df(uS) , (3.15) 

and during the time interval [0,Ti), ut follows the evolution given by the system: dy/dt = F(y). 
At time T\, the momentum is flipped, and gives a new configuration. Then, the system 
starts again with the deterministic evolution up to the time of the next flip, and so on. Let 
£ := (i±, Ti), . . . , (ik,Tj-), . . . be the sequence of sites and ordered times for which we have a flip, and 
let us denote its law by P. Conditionally to £, the evolution is deterministic, and the state of the 
process £jf is given by 

V t G [T k , T k+1 ), w£ = e^- T ^ A o F ik o e (^-^-iM o F ik _, o ■ ■ ■ o e T ' A u , (3.16) 

where F{ is the map u = (r, p) — > (r, p 4 ). 

If initially the process starts from ujq which is distributed according to a Gaussian measure [Xq , then 
&t is distributed according to a Gaussian measure Then, the density fl^ is given by 

£?(•) = J~4(-)dF(0. (3.17) 



D = diag 
D' = diag 



wo 

(3 (x/N) ' PHx/NY 



+ 



Mx/NY 



(3.14) 



More precisely, the mean vector m\ and the correlation matrix of fi\ can be related to the mean 
vector mo an d the correlation matrix Co of {1q : 

ml = e^ A ■ £ ifc • e^-^-iM • • • • e Tl ^ • m , (3.18) 

and 

Cl = e^ T ^ A • • e^~ T ^ A • • • S n • e Tl ^ • C ■ e~ T ^ A ■ • • • e "(^-^-iM . s *. e -(*-T fc M . (3.19) 

Equations (3.18) and (3.19) also give a graphical representation of the process (mt,Ct)t^o- during 
the interval [0,Ti), mt follows the evolution given by the (vectorial) system 

t = F{y) (3.20) 
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(where F has been previously introduced for the process u> t ). At time Ti, the component mj 1+ 7v 
(which corresponds to the mean of p^) is flipped, and gives a new mean vector. Then, the deter- 
ministic evolution goes on up to the time of the next flip, and so on. 

In the same way, during the interval [0, Ti), Ct follows the evolution given by the (matrix) system: 

dM 

(where A has been previously defined). At time Ti, all the components C^ j and Ca x when j ^ i\ 
and i ^ i\ are flipped and the matrix Ct x becomes • CVi • E| . The generator of this Markov 
process (mt,Ct)t^o is exactly the one defined by (3.6). Consequently, for t ^ 0, the law of the 
random variable (mt,Ct) is 9^ no Co , where 



-AM + MA 



(3.21) 



Ao [0 




and 



Cn 



D N 
On D' 



where < 



D = diag 
D' = diag 



1 



+ 



Ag(x/JV) 
Po(x/N) ' /? 2 (x/iV)' 
1 



'Po(x/N)' 

as it can be deduced from Lemma 1.1. Recall that in this section, /j,q is given by 



Z(A)(-),Ao(0) 



It follows that the density /i^ is equal to 

(■) = / Af(-) rfP(0 = / G m ,c(-) dC,co(^.c) • 



(3.22) 



(3.23) 



(3.24) 



(3.25) 
□ 



Remark. Observe that 



Ik IP* 



Af N = / G m , c (r-x) <,c„KC) = / Px d^ iCo (m,C) = E mo , Co [^(t)] 



Gm,c(px) d0* Co (m,C) = / 7r x d^ 0jCo (m,C) = E m0i c [vr x (t)] , 



(3.26) 



(3.27) 



Lemma 3.2. Let (mt,Ct)t^o i>e the Markov process defined above. As previously done, we introduce 
p(t),ir(t) € R N and U(t),V(t), Z(t) € 2%(R) such that 



m t = (p(t)Mt)) andC t =^ ) ^) 



Then, 



\no,C - o. 5. , V t > 0, • 



7Tj(t) < Vy,y(t) , 

P 2 y(t) < ^(i) • 



(3.28) 



(3.29) 
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Proof. First of all, let us notice that the quantities V yy (t) — Tr y (t) and U y ^ y {t) — p y (t) are the diagonal 
components of the symmetric matrix St '■= rat ^m-t — Ct- From Lemma 3.1, we have 

(3.30) 



S t = / Sf dP(0 



For any sequence of sites and ordered times £ = (i\,Ti), . . . , {ik, T/%), . . . , the symmetric matrix S\ 
is positive because this is the matrix of covariances of uf. It follows that St is positive, and its 
diagonal components are all positive. □ 

Remark. In the case of the pinned chain, the matrix A is slightly different, but all the notations 
and conclusions are still valid. The initial correlation matrix for the pinned model is not more 
diagonal, but has non-trivial values on the upper and lower diagonals. The initial mean vector is 
equal to 



3.2 The Evolution of (m t , Q 



t t^o 



Thanks to the regularity of /3q and Ao, we know that there exists a constant K which does not 
depend on N such that 



' jj E [(^-) 2 (o) + (Vi,j) 2 (o) + 2(^) 2 (o)] < K , 
^E[Mo) + ni(o)K*, 



(3.31) 



Moreover, one can easily show that 



^£[(M fe (o) + (^) fc (o) 



< K k , for all k > 1 . 



o 



(3.32) 



It results that the two first inequalities of (3.31) are actually uniform in t, in the sense that 
1 



N 



m ,C 



i -j 



ElM*) + 'M*)] 



(3.33) 



We are going to see how this last inequality can be used in order to show (1.17). We denote by Uk(t) 
and i>k(t) the two quantities 



u k (t) = E m0iCo 
Vk(t) = E mo ,c 



E u iM 



E v ^ 



(3.34) 
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Let us make the link with (1.17). In view of (3.25), we can write 



p? 


2k 
Py 


J 


[ G m ,c 






~2k~ 

y 


~ J 


[ G m , t c 


' r 2k 

' y 



d9 



t 

mo, Co 



[m,C) 



(3.35) 



(3.36) 



We use the convexity inequality (a + b) 2k ^ 2 2k 1 (a 2k + b 2k ) - which is true for all a, b G R - to get 

' (Py ~ % + 7T yf k de mo,C { m , C ) 

\ Py -7r y ) 2k ] d6 t m0tCo (m,C) + 2 2k - 1 [ < oA (m,C) . (3.37) 



Pf 



We deal with the two terms of the sum, separately. First, observe that Gaussian centered moments 
are easily computable: 



G 



m,C 



(Py-^y) ={V v , y -rf) 



k (2k)\ 



y "yJ kl 2 k ■ 



Then, 



E 

y^ N 



K.„ - d9l 0>Co (m, C) < | 



'v,y "yj kl 2 k 
In the same way, 



E 

yeT]v 



1141 



kl 2 k 



mo, Co 



E *?(*) 



E 



(3.38) 



(3.39) 



(3.40) 



Lemma 3.2 shows that 



E mo,C 


E 






E 








E 






E 




= u k (t) 



(3.41) 



As a result, 



2k 
1 II 



,2k 



(3.42) 
(3.43) 



In a few words, to get (1.17), we need to estimate the two quantities Ufc(i) and Vfc(i), which are 
related to Ct- That is what we do in the next section. 

Remark. In the case of the pinned model, the p x and q x remain centered during the evolution: for 
all t > 0. m t = K 2jv . This simplifies the study since we do not need to center the variables. The 
result is the same: we need to estimate u^it) and 
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3.3 The Correlation Matrix 



Lemma 3.3. For any integer k not equal to 0, there exists a positive constant K which does not 
depend on N and t such that 

{ b (3-44) 

\u k (t)^K k N . 

Proof. First of all, (3.33) shows that, uniformly in t, 



U!(t)^KN \vi(t)^KN 

, , an d \ , s 

u 2 (t) < KN \ v 2 (t) < KN . 



(3.45) 



We observe that 



Uk(t) +v k (t) = E m0jCo 



E c 5w 



ieT]> 



(3.46) 



Thanks to the dynamics description, we know the expression of the correlation matrix: conditionally 
to f, for all t £ [T k ,T k+1 ), 

C*(t) = e^ T ^ A -E ik . e ( T *- T *-iM • • • E fl -e^-Co-e-^-Si* ■ • • e-^-^-O^E^-^M , (3.47) 

Consequently, since Co and C^(t) are similar, we have: 

V fc G N, Tr([C^{t)] k ) = Tr(C fc ) = O(N) . (3.48) 

More precisely, 



- E ^(o) + W» - E ^ + + «)' 

From (3.31) we get Tr(Cg) ^ NK k , where K does not depend on N, £ and t: 



(3.49) 



K : = sup 



1 Ag(u) 



ue[o,i] /5q(«) 



(3.50) 



Now we show that the same inequality holds for ^JCf J fc (i). The matrix C^(t) is symmetric, hence 



diagonalizable, and after denoting its eigenvalues by Ai, ...,\2N, we can write 

2N 



Tr([C«(i)] fe ) = E A *- 



(3.51) 



i=l 



We have now to compare X^i=i ^ w hh Y^=i[^ii\ k (^)- But, if we denote by P the orthogonal matrix 
of the eigenvectors of C^(t), then we get C^(t) = (Pf)* ■ D ■ pf, where D is the diagonal matrix with 



27 



the eigenvalues Ai, A2i\r. For the sake of simplicity, we denote by {Pij)ij the components of Pj*. 
Then, 

Hi]Ht)=^EPt j D j ,iPi,}j ^/',,A 7 /' ; ,j =|^-A,| • (3-52) 

But, ^2jPtjPj,i = 1) since D is an orthogonal matrix. Consequently, we can use the convexity 
inequality, and we obtain 

DC&]*(*) < EEW* < E A i = MCW) < ■ (3-53) 

» » J 3 

Hence, 

+ «*(*) < /" JV-K"* < WIT* • (3.54) 

□ 

Remark. We notice that the same proof works for the pinned case. The only difference is about the 
initial matrix Co, but the smoothness of the profile /3q is still true, and the estimate Tr(Cg ) = O(N) 
is valid. 



3.4 When [iff is a Convex Combination of Gibbs Measures 



As in Theorem 1.2, we now suppose that the initial probability measure [iff is a convex combination 
of Gibbs states defined by 

= f G mo ,c (-) da(m ,C ) . (3.55) 



If initially the process starts from ujq which is distributed according to a Gaussian measure G mo c , 
we know from Lemma 3.1 that Cbt is distributed according to a convex combination of Gaussian 
measures written as 

Consequently, in the case where Hq is given by (3.55), the law of the process Cot is given by 



Af(-) 



G m ,c(-) d^ 0!Co (m,C) \ da(m ,C ) . 



Let us recall that we want to control, for k 1, jlf SxeT Px 



and jit 



E 



(3.57) 
. Following 



the lines of the previous section, we notice that it is sufficient to control two quantities: 

da(m ,C ) , 

da(m ,C ) . 



J E mo,C*o 


E U "M 


< 






E 







(3.58) 
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Lemma 3.3 gives a constant C(Aq, A)) which does not depend on N and t such that 



More precisely, 



E m ,C 


E u U f ) 




ief N 


< 




E mo,C 


E 








) = sup <^ 




UG[0,1] I 



1_ AgCtO 



(3.59) 



In order to keep the same control, we have to suppose that, for all k ^ 1, 

/ [K(m, C)] k da(m,C) < oo, where K(m,C) := sup Q 



(3.60) 



(3.61) 



Finally, let us observe that all estimates are given for fi^ but are still true for the accelerated law 
fif. Indeed, the constants that appear do not depend on N and t. 

A Proof of the Taylor Expansions 

Now we prove Proposition 2.1. For the sake of simplicity, we define 

g x (r,p) := -^(p x+l + p x + ^-r x ) , 
27 z 



/x(r,p) := 



Px+l 

7 



(A.l) 



k(r,p) 9x + K(^) f* = F (t, % ) ■ r x h{v, p) . 

First we will compute the first part that appears in the integral N 2 (0^) 1 C* N <f>^ , then we will 
compute the second part —dt4>^ / <j)f x . 



A.l First Term: the Adjoint Operator 
Lemma A.l. 



E # 



xGTjv 

lN 



N 



p x +ir x + 



27 



A? 



2V 



r x +i 

7 



+ |i E [^(^) + ^)] + o(l). 



(A.2) 
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Proof. First, remind that the expression of <j>^ is given by 



Sf(r,p) 



z(xt(-)) 



[ exp I e x (-(3, 



N 



N 



1 



1 -^[^ r x + -F(t,- -r x h(r,P 



x 



N \ N 



By definition, 



a iN iN 



We write down the two conservation laws: 



Me x ) = j x+1 - j% where j% := p x r x _ x , 
A(r x ) = f x+l - f x where f x := p x . 



Hence, 



E K 1 - a (I)) voi). - a, (I) vox),] + % E 

xGTjv 



(A.3) 
(A.4) 



(A.5) 
(A.6) 



(A.7) 



xGTjv 

where V(/) x = f x+ i - f x . 

We are interesting in the first two terms in the sum, and we compute a discrete summation by part. 
Indeed, 



j/GTjv J/GTjv 



We obtain the following terms: 
' x + I s 



A, 



N 
x + l 
N 



N 



, / x \ 1 



/M-) = #(-)-+/3n-)^+o(^3 



iV 



at; tv 



, / x\ 1 



i\ 1 



' Vat; n \n) n 2 



nj n 2 

X \ 1 



AT3 



(A.8) 

(A.9) 
(A.10) 



First of all, we look at the term obtained in the sum with O [N 3 ) . We want to prove 

W E V^xr x O (±^\ f t N d^ < C H N (t) + o(N) . (A.11) 

xGTjv 

We use the entropy inequality Let e : N — > M. be a bounded function. We get 
jj I E P*+i r * W^S) < ^ + ^ log / exp ( A? E^* e (W) <^?o • (A-12) 



But, let us recall the inequality p x+ \r x ^ {p x+ \ + r z)/2 an d for large enough, we have 



Xt(-) 



( %pMn) 



2ir 2N 



J JV->oo 



Np - 2ae{N) 



(A.13) 
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We obtain a similar estimate for [exp (aN 1 r x e(N)^ . 

Therefore, we have showed 



^ / E <N)ffd^ fi < ^ + O(l) . 



(A.14) 



Hence, 

<4E [* (£) + a : (I) *«] + % E [* (£) p « r - + A " (I) fc 



+1 



+ ^E^>+°(^) 



(A.15) 



xGTat 



Moreover, we can compute two equations which are called "fluctuation-dissipation equations". In 
other words, we decompose the current of energy and the current of deformation as the sum of a 
discrete gradient and a dissipative term: 



p x+ i = V 
p x+1 r x = V 



+ £*(/*) 



pi + r x r x _i 
2 7 



(A.16) 
(A.17) 



We use the two equations (A.16) and (A.17), and we obtain 



V 



Pi + r^r^-i 
2 7 



+ C*(g x 



V 



+ C*{f x 



+ 



^ E (J?) ^ + *.* (|) w] + ^ £ + » 

We sum again by part, on the two terms with a gradient, and we obtain as before 

kN 



(A.18) 



x£T N 



P 2 X+1 + r x r x+ i 
2 7 



+ p x +ir x 



r x +i 

7 



+ Px+1 



We get the result. 
Lemma A. 2. 



(A.19) 
□ 



S4> 



N 9t 



N 



E s ^ + §p E ( E w) - ) + $ < N ) 



(A.20) 



where $ [N 2 e(N)} = o{N) . 
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Proof. Thanks to the exponential term, we have 



1 J2 5 x ( p y)-5 x (p) 



(A.21) 



The main idea consists in noting that e x — 1 = x + x 2 /2 + o(x 2 ). We are going to give a rigorous proof 
of this estimate in our context thanks to the hypothesis on the energy moments. More precisely, in 
view of (1.42) and Lemma 2.1, we want to prove that 



fcj N k 

y£T N k^3 



o(N), where F y = Y] (4(p' y ) - 4(p)) • 



(A.22) 



Let us compute F y . We notice that in the following expression, 



(A.23) 



the only terms which are changing when we flip p into p y are 
• the term when x = y, and the difference is 



r yPy of fV 



-R' — 

7 MiY 



(A.24) 



the term when x = y — 1, and the difference is 



r v-tPy at ( V~ l \ , x ,( V- 1 \ 2 Py 



(A.25) 



In other words, we have to show that 



Nf4 



N 



I J? \k 



iV->-00 



■> . 



(A.26) 



with 



^ C \r y p y \ + Ci |r y „ip y | + C 2 \p y \ 
^ K (l + e y + e y -i) , 



where K is a constant which does not depend on N and t. 
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First of all, we introduce the space A y = {e y ^ 1, e y -\ ^ 1}. 

V- (e y + ey-! + l) k K k l {ey ^ ey _^ 1} 



fc>3 



fc! N k 



(3K) k 
fc! N k 

(3K) k _ 

fc! N k N- 



"EE 



« 2 £ 



-+ (A.27) 



Since we have (e y + e y -i) k t A c < (2e y + e y -i) k , we deduce (e y + e y -\) k \ A c ^ Cq e k + C\^ y _ x . 
Consequently, 



JV 



£ — 



K,\ k K k 1 



AC" 



fcS;3 



fc! iV fe 



A: twA: 
fc! iV fc 

fc>3 



+ N E a 



E 

fc^3 



Now we deal with fif [£ fc>3 e k /(k\N k )} . Remind that e k < 2 (pg* + rf ). 

We are reduced to prove that 

n n „A AT \ ,,A r 



fc TWfc 



fc! iV fc 



. (A.28) 



^ E A 



E 

fc^3 



k\ N k 



and N ^ fjq 

yeT N 



E 

fc^3 



fc! N k 



> . 



Af-s-oo 



(A.29) 



We can flip the summations thanks to Fubini theorem. From the hypothesis on the moments bounds 
we get 



yer N 



N 



:2k 



k^3 



fc! N k 



^ 2 E 



fc>3 



(C fc) fc _ 
fc! iV* jv- 



. 



(A.30) 



This last limit is deduced from the property of the series S(x) := X)fc>3 k k x ^ 2 1 (^0 ■ It is a power 
series which has a strictly positive radius and is continuous at 0. Then, 



N ^ { -^M = c 2 s(^ 



fc>3 



fc! N k 



N J N- 



+ . 



(A.31) 



The same happens for the second sum. It follows that 



^Ea 



E 

fc>3 



pk 

y 

fc! A^ fc 



-» . 



(A.32) 
□ 



After adding the two terms and get some simplifications, we obtain this following final result. 
Proposition A.l. 



— N 2 C* rh N 
1\ L N (p t 



e{-^(«4 



pl+i + r x+ ir x 
27 



+ p x+ ir x 



( x 



r x +i 

7 



+ Px+l 



t, + K ^— ^ + 2<9 g A ft, ^1 



+ o(N) . 
(A.33) 
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Proof. There are simplifications when we write (— A + 7>S)(</>£ ). Actually, 

£ {-A(5 X ) + 7 S(5 X ) - C*(6 X )} = . 



ieT f 



The result follows. 



(A.34) 
□ 



A. 2 Second Term: Logarithmic Derivative 

First, we notice that dt<j>^ I4>t = ^*{l°g(0t )}• Moreover, 



log(^) = c + £ „ (-ft (i) + 1) - A, (±) r. - A (i) + P, + f .) 



iV 



iV7 2 7 7V 



+a :(I)^- 1o «[ z < a <-»- a '"»i 



(A.35) 



We need to estimate the partition function Z(fit(-), At(-)). More precisely, we compare this new 
partition function to the exact partition function 



(A.36) 



We prove the following lemma. 
Lemma A. 3. 



cVogZ(A(-),A t (-))-d 4 logZ(A(-),A t (-)) =0(1) whenN^oo. 



(A.37) 



Proof. First of all, remind that the exact expression of Zt '■= Z(j3t(-), At(-)) can be written as 



Z t 



iV 



1 q> ( x \ r ' x ( , ,7 \ 1 W / E \ Px+1 

"n* iivJ ^ + ^ + 2 r -J " a7 a * IaTJ — 



<ipdr 



exp |1||6 £ || 2 | jf ^ exp |-i(A - b t ,C t (X - b t ))jdX = exp !^\\b t || 2 | (2^1 det(C* 



1 1/2 



(A.38) 



where 6t is a vector and Cj is a symmetric positive matrix. 
More precisely, one can see that 



(A.39) 



where ht is a function that can be easily expressed with \t,f3t,\' t and j3[. Then, ht is smooth. 
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Moreover, Ct can be written as Ct = D t + N 1 H t with D t a diagonal matrix and H t a symmetric 
matrix which has at most three non-zero components on each row and each column. More precisely, 



Di 



Hi 



(■. 

\(0) 

V(o) 



\(0) 



(0) 



(oA 



Pt(x/N) 

(•: 
\(0) 



(0) 



(A.40) 



-(2 7 )-^' t (x/N) 
-(2 7 )-^' t (x/N) 



■J 

-(2 7 )^^(x/N) (0)\\ 
-(2 7 )" 1 ^(x/iV) ■-. 

•■/ 

(0) 



(A.41) 



Now we write 



dt log Z t = \ £ d t (| (£) ) + ift log det(Q) + 1 £ ^ 

a tlog i t 4 Ea (|(|)) 



+ -a t logdet(A) • 



(A.42) 
(A.43) 



But, \N~ 2 J2x d th t (x/N)\ = 0(1) since h t is smooth. 

It remains to show that the following quantity is bounded above by a constant that does not depend 
on N: 



dt log 



det C t 
det A 



logdetfl+^A" 1 ^ 



^{det^ + A" 1 ^/^)} 



det(/ + D^Ht/N) 



(A.44) 



We denote by iff the matrix D~j~ l Ht, which also has at most three non-zero components on each 
row and each column, and by K[ the derivative of Kt with respect to t. We notice that for N large 
enough, the matrix I + Kt/N is invertible, and we have 



d t log 



detCt 
det D t 



Tr( *com(J + K t /N) ■ (I + K'JN)) 



det(/ + JTt/W) 



Tr 



I+—K t 

N 



il + K't) 



where com(^4) is the comatrix of A. 
Now we deal with (I + Kt/N)- 1 : 



(A.45) 



(A.46) 



But, the component (i,j) of K* can be written as i k ai M a u,«2 ' ' ' a ifc,j wnere a i,j are the 

components of Kf. We know that there are at most three non-zero components on each row and 
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each column, and that they are all bounded by a constant C that does not depend on N (since j3 t 
and At are smooth). Then, it implies that |Tr(ET t fc )| N3 k C. 

It follows that 



Tr 



< 1 + wi+ c E^ = °( 1 



k^2 



(A.47) 



because Tr(ATj) = O(l) (we can compute it and again use the smoothness of the profiles). 
In the same way, we show that 



Tr 



K' t [l+-K t 



0(1) . 



It ends the proof. 



(A.48) 
□ 



We deduce from the previous result that 



dt log [Z (p t (-),A f (■))]= 



d t Pt(x/N) 
' Pt(x/N) 



+ d t \t(x/N) 



X t (x/N) dMx/N) X 2 t (x/N) 
Pt(x/N) 2 ${x/N) 



+ 0(1) . 
(A.49) 



Consequently, we get the following statement. 
Proposition A. 2. 



7 A 



x\ dtP(t,x/N) 



NJ P(t,x/N) 



dtX u x/N) MhlM + W,x/N) X 2 (t,x/N) 



P(t,x/N) 



p 2 {t,x/NY 



(A.50) 



ei-e U 



A" 



dtP(t,-) + 



r x -r[t, 



N 



dtX[t,-)+0(l) . (A.51) 



A. 3 Ending Proof 



We are now able to prove the Taylor expansion. According to the results of the two previous parts, 
we have 



^N 2 C N ^ - d t {log(^)} = £ \ -d\p (t, ^ 



a;£T» 



2 7 



r x . 
VVx 

7 



+ 



[(r x + r x -i)d q p(t,^) 



r r — r ( t, — 

' N 



d t x(t,^)} + o(N) 
(A.52) 



3G 



Using the notations introduced in Section 2, it becomes: 

-Liv^f _ 9t{log{(j) N )} = E r_ i q p ^ |) ji _ i d 2 X ^ |) j; 



+ 



47 
1 

+ - 

7 



2 4 + ~d q P (t, |p) 5 9 A (t, |) 4 



+ 



x 



e x -e[t, 



N 



1 



x 



9tP(t, w ) 



r x -rlt 



N 



dtXlt, 



N 



)} + o(N) 
(A.53) 



We denote by Hf. the function defined as follows: 



Hk[ri[t,- 



J 



(A.54) 



The explicit formulations for are given by Proposition 2.1. The sum 



2^1^) 



.r 



1 

—i 

1 



A/ 



cl 2 A It, — )H 2 [v [t, 



N 



1 

1 

+- 

7 



iV 



1 



N 



N 



x 



N 



2 / / X 



(A.55) 



is of order o(iV) (thanks to the regularity of the functions e, r, (3, A), so that we can introduce it in 
the right member of the equality (A. 52). 

Then, we obtain after computations 
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and 
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drH, - dr H 2 + M- «9 r tf 3 + 



<9 g /?d g A fl ^ , [3,A]' 



47 7 
Indeed, these two quantities are respectively equal to 
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(A.58) 



and 



ffi r , d 2 q X [d q (3] 2 
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r (2e - 3r 2 ) - ^ (2e _ 3^) + r M 



7 



This concludes the proof and gives Proposition 2.1. 



(A.59) 
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B Proof of the One-block Estimate 



We just give a sketch of the proof, which is done in [4], Section 3.4. First, we define the space time 
average of distribution: 

(B.l) 



i=i j[) 



and f* its projection on {(r^Pi) G M 2 ( fc+1 ) ; i g A k ■= {-[jfe/2] - 1, ... , [k/2] + 1}}. 

We also denote dv N = f N Y\ i€ j N dridpi and dv^ = f^ N \\ i( z.f N dridpi the corresponding probability 
measures on R 2N and R 2 ( fc+1 ) . 

Observe first that (2.22) can be rewritten as 



t lim sup lim sup lim sup 

M— >oa £— >-oo N~>oo 



iGA £ (0) 



dv N = , 



(B.2) 



because 

1 t-l 1 P 

fc=0 ^ j=l 

We can prove the first following lemma. 
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(B.3) 



Lemma B.l. For each fixed k, the sequence of probability measures (i / L^)jv>fe is tight. 

For any k let be a limit point of the sequence {v^)n^i- The sequence of probability measures 
(ffc)fcj>i forms a consistent family and by Kolmogorov's theorem there exists a unique probability 
measure v on (M x R) z such that the restriction of v on {(ri,pi) € R 2 ( fe+1 ) ; i £ A&} is u^. One has 
easily that v is invariant by translations. 

Lemma B.2. For any bounded smooth local function F(r,p), we have J CFdv = 0. 



Then, v is a convex combination of grand canonical Gibbs measures p x 
with p a probability measure such that J dp(x)p x [ej] ^ Cq for any j € Z. 

Hence, it results that 
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3S 



where the last equality is a consequence of the dominated convergence theorem. Since [i x is ergodic 
with respect to {t x ; x € Z}, the last term is equal to 

limsup / dp(x) \n x [J ,m] - #(Mx[ ? ?o,m])| • (B.5) 

As M -> oo, h x [Jq,m\ converges to /U X [J ] = # (ju x [£oD and /x x [£ 0) m] to /x x [£ ]- 

By Fatou's lemma, the limit in M is equal to and this concludes the proof of the one-block lemma. 
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